Wentzel-Bardeen singularity and phase diagram for interacting 
electrons coupled to acoustic phonons in one dimension 
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We consider strongly correlated electrons coupled to low energy acoustic phonons in one dimension. 
Using a Luttinger liquid description we calculate the exponents of various response functions and 
discuss their striking sensitivity to the Wentzel-Bardeen singularity induced by the presence of 
phonons. For the Hubbard model plus phonons the equivalent of a phase diagram is established. By 
increasing the filling factor towards half filling the Wentzel-Bardeen singularity is approached. This 
in turn triggers a simple and efficient mechanism to suppress antiferromagnetic fluctuations and to 
drive the system via a normal metallic state towards a superconducting phase. 
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Recently there has been a renewed interest in one- 
dimensional (ID) interacting fermion systems, both in 
the context of high temperature superconductivity [Q, 
and in the context of mesoscopic systems (^-|5). ID mod- 
els are easier to handle than their 2D and 3D counter- 
parts, and can yield useful information for the latter sys- 
tems. It has long been recognized that in ID elec- 
tron systems the low-energy long wave-length excitations 
around the Fermi surface have bosonic properties. The 
electron system is then represented as an harmonic fluid 
[@-||, called Luttinger hquid [0, where the Fermi oper- 
ators are expressed in terms of scalar Bose fields. 

In this letter, we study the properties of a ID system 
of correlated electrons, interacting with a short range re- 
pulsive potential, and coupled to small-momentum trans- 
fer bosons with acoustic dispersion. We concentrate 
first on acoustic phonons, and later on discuss an ap- 
plication to electron waveguides. While the low dimen- 
sionality and the short range forces rule out the possi- 
bility of long range order at any finite temperature T, 
one can study which type of ordering fiuctuations domi- 
nates the phase at vanishing T - charge/spin density 
waves (CDW/SDW) or singlet/triplet superconducting 
(SS/TS) fluctuations - and plot the equivalent of a phase 
diagram. 

ID electron systems coupled to small-momentum 
phonons were considered more than four decades ago by 
Wentzel |l^ and by Bardeen jl2j as candidates for the 
theory of conventional superconductivity. It was pointed 
out that for a critical phonon coupling constant, the sys- 
tem becomes unstable, and the specific heat diverges as 
one approaches this critical point JT^ : the retarded in- 
teraction mediated by the phonons induces a collapse of 
the system. Below, we refer to this singular point as 
the Wentzel-Bardeen (WB) singularity. Recently, an ex- 
actly solvable model similar to that of Ref . 1 1^ was stud- 
ied urn , with the conclusion that these small- momentum 



phonons can be ignored in typical metals since their effect 
is of the order c^/vp <C 1 {vp is the Fermi velocity and 
c the sound velocity) . Consequently, the interest jlj-|l^] 
shifted towards larger-momentum {2kp) processes, where 
an electron is scattered by a phonon across the Fermi sur- 
face. 

However, we shall show that the small-momentum 
phonons do play an important role in strongly corre- 
lated systems, as it is the velocity Up of particle-hole 
excitations which determines the "adiabatic parameter" , 
(? jv^p, and not the Fermi velocity vp. In particular, we 
argue in the following that: 1) the coupling to small- 
momentum transfer modes (such as acoustic phonons) is 
a non-perturbative effect as one reaches the WB singular- 
ity. 2) For the Hubbard model, the WB singularity can 
be reached for arbitrary electron-phonon coupling con- 
stant by increasing the filling factor towards half filling. 
3) Near this WB singularity, CDW and SDW fluctuations 
are suppressed, and the system is pushed into a metallic 
and finally superconducting phase by slightly increasing 
the filling factor. Migdal's theorem does therefore not 
apply in the presence of strong correlations. Finally, it 
is worth pointing out that this high sensitivity towards 
doping is reminiscent of high temperature superconduc- 
tivity materials. 

We use the Luttinger liquid description generalized to 
electrons with spinjq], later on to be further specified to 
a Hubbard model p9| |. Our starting point is the Hamil- 
tonian, H = He + Hph + iJe-ph, where 



(0.1) 



describes the free phonons, with d the displacement field, 
lid its canonical conjugate and C, the mass density. The 
electron-phonon coupling is 120] : 



1 



H, 



el-ph — g 



(0.2) 



where g is the coupling constant, and tps = 
^ikpx^^^ _l_ ^-ikpx^^_ ^Yie electron field operator. In 

terms of boson fields, ■0s±(x) c>c exp[iY^(±(^s(x) — 
/'^ ns(x', r)dx'))] (here x = {x,t)) representing right 
and left moving electrons. The electron density in Eq. 
( |0^ reduces to Y^s'^l'^^ = \f^dx^p., with ipp^„ = 
((/?! ± (/?|)/V2 the charge and spin fields (similarly for 
the canonical conjugate fields Ilp cr). We have neglected 
the fast oscillating terms in the density which give rise 
to 2fci? phonon induced backscattering. The electronic 



Hamiltonian is 



H„ + H„, with 



dx 



(0.3) 



with V = p,a. The phonons couple only to the charge 
degrees of freedom, and this property is preserved in the 
propagators calculated below. In Eq. (^), the elec- 
tron interaction is included in the parameters Kp, Kg-, 
and the charge (spin) velocities Up (ua)- This model be- 
longs to the same universality class as the multicompo- 
ncnt Tomonaga-Luttinger model jl^ . 

To calculate Green functions, we use a Lagrangian for- 
mulation, where the partition function is represented as a 
functional integral. As we will be calculating Green func- 
tions which involve fermion operators only, the phonon 
degrees of freedom are integrated out right away. This 
leaves us with an effective action for the charge de- 
grees of freedom, which reads in Fourier representation: 

s, 

inverse propagator 



p - (2/3L)-iEk^p(k)l¥'p(k)P (k = {k,u;)), with the 



^p(k) = 



KpUp 



2 I 2i 2 

oj + Upk 



(0.4) 



where b = g^jKpUp/C,. The first two terms of Eq. (0.4) 
represent the contribution of the free charge field, and the 
retarded, attractive coupling associated with the phonons 
appears in the third term. At the WB point the charge 
density propagator, fc^/I?p(k), becomes proportional to 
w"^ and signals an instability towards long wave-length 
density fluctuations. 

We first calculate the single-particle Green function: 



Gs{x,t) = -< rVa(a;,T)V'l(0,0) > 



(0.5) 



where T is the imaginary time ordering operator. This 
quantity has been calculated before pl[ | (by diagram- 
matic methods) for spinless fermions coupled to acoustic 
phonons, but in the perturbative limit /vp <C 1, where 
the WB singularity is absent. Using the decomposition 
of the Fermi operators into right and left moving compo- 
nents, and after normal ordering |^,^, the calculation of 
the Green function is reduced to an evaluation of Gaus- 
sian integrals. The result is (a is the lattice constant): 



Gsix,T) 



2n 



a 



a 



xn 

X e 

Q = ± 



X + lUo-T 

a 



KpUpF^/Avp+Cp/AKp 

X + ivpT 

, ■ I I \ sqn(T)/2 
sax + tUcr\T\ \ " ^ " 



n 



ax + iVj\t\ 
i\x + iv-yTl 



i\x + iu„t\ 

sgn(T)F^/2 



7=± ^ ' ' 

xna{x,T) , 
with the normal ordering contribution: 
n„(a;,r) =sgn(T) e-""^s"M[«-(^'" 



(0.6) 



),V„(0,0)]/2 



(0.7) 



and thus = 1. If x/L 7^ and t arbitrary, ^^(x, t) = 
sgn(T) exp(iQ;7ra;sgn(r)/L), and if r = 0^ and x arbi- 
trary, na{x,0^) = ±ia\x + iaa\/{x + iaa). Near x = 
one has to add a to v\t\ in Eq. (^). We introduced 

the velocities, v\ = [v? + =F 



^(u2_c2)2+462]/2. 



for the electronic collective excitations. These veloci- 



ties characterize the roots of the propagator (0^), and 
arise from the hybridization of the charge and phonon 
excitations. Note that it is the charge velocity Up (not 
vp) which appears here. The exponents are given by 



C± = Up{v\ - + b 



V<)/[(4 



i;2:)i;±], and F± 



± 



At large distances, the Green function decays as a 
power law, Gs{x,0) oc with S ~ K„/A + 

l/(4if„) - 1 + B/{AKp) + AKp/A. The electron-phonon 
parameters are defined hy A = Up{l + /v+V-)/{v+ + 
v^) > 1, and B = Up{l + v+v^/ul)/{v+ + w_) < 1. 
These parameters play a crucial role in the discussion of 
the ordering fluctuations below. In the limit g ^ 0, 
A = B — 1. Let us now consider the case g 7^ 0, 
where A > 1 and B < 1 |Q: As Up/Kp approaches 
the critical value u*p/Kp = g"^ /{C,c^) from above, 



tends to zero, and v'^_ to 



c^. As a result, the ex- 



ponent A diverges and B decreases to the finite value 
B* = u*pl^uf + c2 < 1. For Up/Kp < u*p/K*p the ve- 
locity becomes complex and the model becomes un- 
physical. Thus we must require that Up/Kp > u*/Kp, 
or equivalently that b/ (cu p) < 1, the equality sign de- 
fines the WB singularity ||23|| . This singular behavior 
is a non-perturbative effect of the electron-phonon cou- 
pling |24[| an d originates from the instability of the prop- 
agator (Q. Because J > 0, if g 7^ [§2), the Green 
function decays faster than l/x, the free electron result. 
From ( p.6D , we find the momentum distribution function 
near kp, N{k) ~ N{kp) — KSgn(/c — kp)\k — kp\^ , with 
K some constant of order one, and N{kp) = r(l/2 -|- 
i5/2)/[2-y7rr(l + (5/2)]. This continuous but non-analytic 
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behavior at the Fermi surface is characteristic for Lut- 
tinger hquids . The presence of phonons induces this 
Luttinger hquid behavior even without electron-electron 
interaction {Kp = — 1). 

We now turn to the discussion of the ordering fluctu- 
ations characterized by two-particle response functions. 
The definitions for N{x), x(x), As(x) and A((x) which 
describe CDW, SDW, singlet (SS) and triplet (TS) su- 
perconducting fluctuations are given in (|,^. The calcu- 
lation of these Green functions will be given elsewhere 
p2[ . Here, we give the final result in Matsubara repre- 
sentation: 



(0.8a) 



N{x,t)(x Y[ \x + ivpT\-'^''^^^'''''''\x + iu„T\-'^'' 
X(x,r)(x W \x + ivpT\-''''^''^^'''''\x + iu^T\-^' 

p=± 

(0.8b) 

As(x,r)cx Jl \x + ivpT\-^^''^''\x^iu^T\-^'' (0.8c) 



At{x,T)(X Y[ \x + iVi3T\ ^f^^''\x + iUaT\ 



(0.8d) 



/3=± 



The signature for a particular ordering fluctuation to be 
present is given by the divergence of the Fourier trans- 
form of the corresponding response function at low fre- 
quency and small momentum relative \o q — lkp (g = 0) 
for TV and x (^s *)• We thus obtain the following criteria 



KpA 



K^<2 



KpA + l/K^ < 2 
B/Kp + K,<2 
B/Kp + < 2 



(CDW) 
(SDW) 
(SS) 
(TS) 



(0.9a) 
(0.9b) 
(0.9c) 
(0.9d) 



Since A > 1 and S < 1 for g 7^ 0, we see that the Cooper 
instability is always present for non-interacting electrons 
coupled to phonons. The exponents can in principle be 
recovered in | p8[ , which treats formally an arbitrary num- 
ber of coupled charge fields: here we have the advantage 
of having explicit analytic expressions. 

We now specialize to a Hubbard model for the corre- 
lated electrons. In this case |19 , Ka- = 1, and the CDW 
and the SDW (the SS and the TS) response functions 
have apparently the same exponents (except at half fill- 
ing). However, logarithmic corrections due to marginally 
irrelevant operators in the spin channel favor SDW (TS) 
over CDW (SS) fluctuations . The remaining param- 
eters Up and Kp of the Luttinger liquid Hamiltonian are 
obtained for arbitrary on-site repulsion U and flUing fac- 
tor n by solving two integral equations [p^|j2^ ] describ- 
ing the ground state properties and the spectrum of 
charge excitations p7|] . 

Our results are plotted in Fig. |l| a) and b): for fixed 
U/t (Fig. H a), where t is the bandwidth, we determine 



which fluctuations dominate as a function of n and an 
effective electron-phonon parameter b/{upc): for conve- 
nience we consider in these plots n (U) and h/(upc) as 
independent parameters. For small b/{upc), SDW (i.e. 
antifcrromagnetic) ordering fluctuations dominate for ar- 
bitrary flUing factor. As this parameter is increased, we 
reach (away from half flUing) a region for which no re- 
sponse function diverges: we refer to this "phase" as the 
metallic region. At low U and b/{upc) ^ 1 analytical 
results confirm the existence of this intermediate phase 
p2[ . By further increasing b/{upc), the region where su- 
perconducting fiuctuations dominate is reached. On the 
other hand, near half filling, the correlation effects sup- 
press the superconducting phase. (For larger values of 
U, the superconducting region shrinks towards the re- 
gion where b/{upc) = 1, and the SDW region grows as 
expected [^.) Next, we choose quarter filling [n — 1/2), 
and plot in Fig. |l| b) the phase diagram as a function of \J 
and the phonon coupling. At low J7, Kp k, 1, and the sys- 
tem is superconducting because of the Cooper instability. 
As U is further increased, one crosses the metallic and 
the SDW (CDW) phase, because the phonon-mediated 
attractive interaction is overcome by the instantaneous 
repulsion between electrons. 

In both Figs. ^ a) and b), the upper line b/{upc) = 1 
corresponds to the WB singularity. We now discuss the 
relevance of this singularity for the Hubbard model. Plot- 
ting Up/ Kp for several values of [/ as a function of the fill- 
ing factor (Fig. ^) , we notice the well-known fact ||l9| that 
Up/ Kp vanishes as one approaches half filling (and, of 
course, zero filling). Thus, the WB singularity at Up/K* 
can be reached for arbitrary values of the coupling con- 
stant g (^ 0). Near half filling, as one approaches u*p/K* 
from lower filling factors, the divergence in A triggers a 
dramatic suppression of the SDW (CDW) fluctuations, 
and the system is driven into the metallic phase [ p8| . 
Moreover, since B ^ B* as Up ^ u*, and B* cx g^, the 
SS/TS condition B < Kp can be met near half filling for 
sufficiently small 5 7^ (since 1/2 < Kp < 1). Hence, 
the system can finally be driven into the superconduct- 
ing phase by approaching half filling. In summary, the 
conjunction of the coupling to small momentum acous- 
tic modes with the reduction of the charge velocity near 
half filling provides a very efficient mechanism to sup- 
press SDW (CDW) order, and to drive the system into a 
metallic and finally superconducting phase. This mecha- 
nism need not be limited to the ID Hubbard model, but 
could occur in other situations where strong correlations 
play an important role. 

We conclude with an application to mesosco pic q uasi- 
ID wires. The phonon field in Eqs.(3.1) and (0.2) now 
represents the charge excitations of another transverse 
state. This model was used in Q to describe the trans- 
port properties of a multimode wire [p9| . The Fermi 
velocities associated with the two modes are compara- 
ble and therefore the retarded attractive interaction is 
much stronger than for phonons. Beyond a critical cou- 
pling, superconducting fiuctuations are induced first in 
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the branch with the lower Fermi velocity, and then in 
both modes. This could be probed e.g. by studying 
the periodicity of the persistent current of a two channel 
mesoscopic ring in a heterostructure |]30| , as a function of 
electron density. The flux causes a twist in the boundary 
conditions, but the bulk properties are unaffected. The 
mutual inductance of the two modes is negligible because 
of the small ratio avp /ci {a ~ 137 and q is the velocity 
of light Q). In the metallic regime, the power spectrum 
has a peak associated with the flux quantum 0o = he/ e, 
whereas in the superconducting case, this peak should 
give place to a peak at the superconducting flux quantum 
(^0/2 

ACKNOWLEDGMENTS 

We thank A. J. Leggett, CP. Enz, and S. Trugman 
for useful discussions. The work of D.L. is supported by 
NSERC of Canada. 



[20] A. Fetter and D. Walecka, Many Particle Physics 
(McGraw-Hill, New York, 1969). 

[21] M. Apostol and L Baldea, Phys. Lett. 88A, 73 (1982); 
J. Phys. C 15, 3319 (1982). 

[22] Th. Martin and D. Loss, preprint (1993). 

[23] This singularity was left unnoticed in Ref. [0 (only for- 
ward scattering is includede there), and the exponents 
given there seem to be incorrect. 

[24] Note also that the limits g —* Q and Up — » u* do not 
commute. 

[25] T. Giamarchi and H.J. Schulz, Phys. Rev. B 39, 4620 
(1989). 

[26] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 
(1968). 

[27] C. F. Coll, Phys. Rev. B 9,2150 (1974). 

[28] Note that the model becomes unphysical when we go 
below the WB singularity u*p/K* in Fig. 2 by increasing 
(lowering) the filling factor towards half (zero) filling. 

[29] For further properties of coupled chains see also M. Fab- 
rizio, A. Parola and E. Tosatti, Phys. Rev. B 46, 3156 
(1992); A.M. Finkelstein and A.I. Larkin, Phys. Rev. 
B47, 10461 (1993). 

[30] D. Mailly, C. Chapelier, and A. Benoit, Phys. Rev. Lett. 
70, 2020 (1993). 



* Permanent address as of October 1994: Institut Laue 
Langevin, B.P. 156, 38042 Grenoble, France. 

[1] P. W. Anderson, Phys. Rev. Lett. 67, 3844 (1991). 

[2] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 
1220 (1992). 

[3] D. Loss, Phys. Rev. Lett. 69, 4630 (1992). 

[4] D. Loss and Th. Martin, Phys. Rev. B. 47, 4619 (1993). 

[5] K. A. Matveev and L. L Glazman, Phys. Rev. Lett. 70, 
990 (1993). 

[6] D. C. Mattis and E. H. Lieb, J. Math. Phys. 6, 304 
(1965). 

[7] A. Luther and L Peschel, Phys. Rev. B 9, 2911 (1974). 
[8] V. J. Emery, in Highly Conducting One Dimensional 
Solids, J. T. Devreese, R. P. Evrard and V. E. van Doren 
eds. (Plenum, New York, 1979), p. 327. 
[9] J. Solyom, Adv. Phys. 28, 201 (1979). 
[10] F. D. M. Haldane, J. Phys. C 14, 2585 (1981); Phys. Rev. 

Lett. 47, 1840 (1981). 
[11] G. Wentzel, Phys. Rev. 83, 168 (1951). 
[12] J. Bardeen, Rev. Mod. Phys. 23, 261 (1951). 
[13] S. Engelsberg and B. B. Varga Phys.Rev. 136, A1582 
(1964). 

[14] J. Voit and H. J. Schulz, Mol. Cryst. Liq. Cryst. 119, 
449 (1985). 

[15] J. Voit and H. J. Schulz, Phys. Rev. B 34, 7429 (1986); 

36, 968 (1987); 37, 10068 (1988). 
[16] G. T. Zimanyi, S. A. Kivelson and A. Luther, Phys. Rev. 

Lett. 60, 2089 (1988). 
[17] J. Voit, Phys. Rev. Lett. 64, 323 (1990). 
[18] K. Penc and J. Solyom, Phys. Rev. B 47, 6273 (1993). 
[19] H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990); Int. J. of 

Mod. Phys. B 5, 57 (1991). 



FIG. 1. a) Phase diagram of the Hubbard model coupled to 
phonons, as a function of filling factor and phonon coupling 
constant b/cup, for the choices U/t — 0.3 and c — at (this 
corresponds e.g. to (c/up)^ ~ 1/2 at quarter filling (n — 1/2). 
SDW, M, and SC label the SDW, intermediate, and (triplet) 
superconducting regions, b) Phase diagram of the Hubbard 
model coupled to phonons at quarter filling (n = 1/2), as a 
function of U/t and b/cup. 

FIG. 2. Plot of Up/Kp as a function of the filling factor 
n. Note the abrupt change for small values of U as one ap- 
proaches half filling. 
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